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Here we report extensive ultrafast time-resolved reflectivity experiments on overdoped 
Bi2Sr2Cai_2:Ya;Cu208+i single crystals (Tc=78 K) aimed to clarify the nature of the 
superconducting-to-normal-state photoinduced phase transition. The experimental data show the 
lack of the quasiparticles decay time divergence at the fluence required to induce this phase tran- 
sition, in contrast to the thermally-driven phase transition observed at Tc and at variance with 
recently reported photoinduced charge-density-wave and spin-density-wave to metal phase transi- 
tions. Our data demonstrate the non-thermal character of the superconducting-to-normal-state 
photoinduced phase transition. The data have been analyzed using an ad-hoc developed time- 
dependent Rothwarf- Taylor model, opening the question on the order of this non-equilibrium phase 
transition. 

PACS numbers: 74.40. -|-k,74. 72. Hs,78.47.J- 
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I. INTRODUCTION 



The possibility of inducing an electronic non-thermal 
phase transition in high-temperature superconductors 
(HTSC), by means of ultra-short laser pulses, will set 
a new path for studying the origin of the supercon- 
ductivity in these materials. In facts, under such non- 
equilibrium conditions, the homogeneous superconduct- 
ing phase becomes unstable as its free energy increases 
during the pulse duration^^"^ , while the superconduct- 
ing order parameter can coexist with the pseudogap or 
the normal-state. For many years the exploration of the 
physics of this process has been a difficult task because 
of experimental limitations, mainly arising from laser in- 
duced heating of the samples''"^. Only recent all-optical 
pump-probe experiments on underdoped and optimally 
doped HTSG^'^'S have achieved control of the impulsive 
vaporization of the superconducting condensate in the 
high-intensity regime. This phenomenon has been ob- 
served as the saturation of the transient reflectivity vari- 
ation (AR/R) signal associated with the superconducting 
phased, in contrast to its linear fluence dependence in 
the low intensity regime^iiSrJl. 

A similar photoinduced phase transition (PIPT) to 
the metallic phase has been recently reported on 
charge-density-wave^^'^^ (CDW) and spin-density wave^^ 
(SDW) compounds. While in such systems the PIPT 
exhibits quasi-thermal character, as it is accompanied 
by the same critical slowing down observed in quasi- 



equilibrium conditions at Tc^"— , the origin of the PIPT 
in HTSC remains unclear. A picture of the photoinduced 
non-equilibrium state is still lacking due to intrinsic diffi- 
culties to disentangle pseudogap and normal-phase sig- 
nals in optimally and underdoped samples^'^, the fin- 
gerprint of the pseudogap phase being the AR/R sign 
change observed above Tc in pump-probe experiments 
when probing at 800 nm wavelengtbi^. 

Here we report pump-probe optical reflectivity mea- 
surements at 800 nm in the high-excitation regime on an 
overdoped Bi2Sr2Cao.92Yo.o8Cu208-i-5 (Y-Bi2212) single 
crystal (Tc=78K). At this doping regime the underly- 
ing normal phase is Fermi-liquid-likei^ and no AR/R sign 
change above Tc is measured, at variance with optimally 
and underdoped samples^^. These characteristics are 
fundamental to quantitatively interpret the data within 
a Rothwarf- Taylor (RT) model with time-dependent pa- 
rameters. 

We show the lack of the quasiparticles decay time in- 
crease at the fluence required to photoinduce the tran- 
sition, at variance with previously reported PIPTs on 
CDW— ii^ and SDW— systems and in contrast to the de- 
cay time divergence observed at Tc on cupratesiSiiiii^. 
This finding demonstrates the non-thermal character of 
the superconducting-to-normal-state PIPT, opening the 
question on the nature of this phase transition. To 
address this question we develop a time-dependent RT 
model. The measured decay dynamics is well reproduced 
by this model, resulting in a non-equilibrium supercon- 



ducting gap of about one half of the equihbrium value 
at the pump fluence threshold for the PIPT. We propose 
two different pictures of the non-thermal phase transition 
occurring in the experiment. 

The present work represents a landmark for the grow- 
ing field of pump-probe techniques, which have been re- 
cently extended to the use of several probes, such as Ra- 
man scattering^l, electron-diffractionS^, angle-resolved 
photo emission--- and broadband optical spectroscopjj^i. 
All these techniques require an intense ultrashort pump 
laser pulse, ranging from 0.5 to several mJ/cm^ (Ref. ,23. 
and [23 respectively), to obtain reliable results. Our re- 
sults tackle the long-standing question on the effect of a 
pump laser pulse at high fluence on the superconducting 
condensate of HTSC. 

In Sec. [Uwe briefly report the experimental procedure. 
In Sec. IIIIAI we study the temperature dependence of the 
signal related to the superconducting phase in the low- 
excitation regime and we show, within the context of the 
RT model-- , that the decay rate is proportional to the 
superconducting order parameter A. In particular, we 
focus on the decay time divergence in the vicinity of Tc, 
where A ^- and thus a vanishing relaxation rate is 
measuredi2iiiiHi2£. 

In Sec. IIIIBI we discuss the dynamics as a function 
of fluence at a fixed temperature (10 K) well below Tc- 
Above a threshold pump fiuence, the refiectivity varia- 
tion deviates from the linear dependence and exhibits a 
saturation in agreement with previous experiments^i^i^. 
This discontinuity is identified with the condensate va- 
porization in the whole probed volume^. We observe the 
absence of a decay time divergence above this thresh- 
old point at variance with the quasi-thermal case (Sec. 
IIII A[) and with experimental observations of PIPT on 
CDWiS'i'' and SDWil compounds. 

To interpret the data in the high-fiuonce regime we take 
into account dynamical variations of the non-equilibrium 
superconducting gap, thus extending the RT model to 
the time-dependent case (Sec. IIIICp . In Sec. IIII Dl we 
present the experimental data within the analytical re- 
sults of the RT model. We obtain a linear decrease of 
the non-equilibrium gap with pump fluence and a finite 
non-equilibirum superconducting gap 2 ps after the pump 
pulse that causes the PIPT. 



E- 2- 



-6- 




-T=95 K 
-T= 91 K 
-T= 86 K 
-T= 82 K 
-T- 78 K 
-T . 74 K 
-T = 70 K 
-T.66K 
"T= 62K 
T=58K 
T= 54K 
T= 50 K 
-T = 46 K 
-T = 42 K 
-T-38K 



~i r 

1 2 3 

delay time (ps) 



FIG. 1: AR/R signal (solid lines) as a function of delay time 
at different temperatures (from 95 K to 38 K) on Y-Bi2212 
overdoped single crystal. The dashed lines are the fits to 
the experimental curves (see the fitting function in the text). 
The black line represents the AR/R trace at Tc. In the 
inset we report the fit parameters of the SCS, i.e., decay time, 
tb (red circles), and amplitude, B (black circles). The error 
bars obtained from the fit procedure are displayed for all the 
data. The high uncertainty on the decay time above Tc led 
us to plot just the error bar and not the data point. The 
black arrow on the temperature axis of the inset shows the 
superconducting transition amplitude of the sample. 



ramie at the same oxygen content. In our experiments 
the 800 nm, 120 fs laser pulses are generated by a cav- 
ity dumped Ti:sapphire oscillator. The use of a tunable 
repetition rate laser source allows us to avoid the ex- 
perimental problem of average heating effects^^. These 
effects prevented earlier observation of the photoinduced 
condensate vaporization^'''. In the time-resolved experi- 
ment the transient reflectivity variation is measured and 
we denote it as AR/R. 



III. RESULTS AND DISCUSSION 



II. EXPERIMENTAL METHODS 



A. Low-fluence results as a function of temperature 



Pump-probe measurements have been performed on 
an overdoped Bi2Sr2Cao.92Yo.o8Cu208-i-5 single crystal 
with Tc=78 K ± 5 K. The Y-substituted Bi2212 single 
crystal was grown in an image furnace by the traveling- 
solvent floating-zone technique with a non-zero Y content 
in order to maximize To^l. The crystal was annealed in 
flowing oxigen in order to increase the hole concentra- 
tion and reach the overdoped side of the phase diagram. 
The sample was subsequently homogenized by further 
annealing in a sealed quartz ampoule, together with ce- 



In the low-intensity regime several common trends 
have been recognized based on the large amount of ex- 
perimental data reported on HTSC: i) the appearence 
below Tc of a AR/R signal proportional to the photoex- 
cited quasi-particle (QP) density^"— i^, ii) an intensity 
and temperature dependent decay time of this supercon- 
ducting component^i^i^"— i^ iii) the equivalence of the 
decay time observed probing at 800 nm and the gap dy- 
namics observed in the Thz spectral region^^. 

In Fig. [T] we plot the raw (AR/R) data of consec- 



utive series of scans taken from 95 K down to 38 K 
at low pump fluence (~4 /xJ/cm^). The positive signal 
above Tf can be reproduced by a single exponential de- 
cay function, yle"*/"^-^, convoluted with the time shape 
of the laser pulse, describing the relaxation of hot elec- 
trons via electron-phonon interaction^Sri^ with a relax- 
ation time TA of about 440 fs. All the traces from room 
temperature down to 95 K did not show any substantial 
difference from the 95 K trace. In agreement with the 
literature^°, no AR/R sign change associated with the 
pseudogap phase is observed when probing at 800 nm 
wavelength on the overdoped sample, at variance with 
optimally and underdoped samplesi^. 

The positive and fast signal survives also below T^, su- 
perimposed to a negative signal that we recognize as the 
superconducting signal (SCS) proportional to the pho- 
toexcited QP density and that we capture with a single 
exponential function. The sum of these two functions, 
^g-ty-TA _j_ Bg-*/'^B ^ jg sufhcient to reproduce all the ex- 
perimental curves as shown in Fig. [1] The relevant fit 
parameters of the SCS, i.e., the amplitude B and the 
decay time tb, are plotted in the inset. 

The decay time of the SCS, tb, is shown, as a func- 
tion of temperature, in the inset of FigH] (red circles). 
Starting from 38 K the decay time decreases with tem- 
perature, reaching a minimum at 62 K of 3.3 ps. Above 
this temperature it rapidly increases reaching the value 
of 4.7 ps at Tc". At this temperature the order parameter 
vanishes (A ^-0), the SCS amplitude, B, is much smaller 
than A (the normal-state signal), and the uncertainty on 
the determination of its decay time strongly increases. 

The increase in the decay time in the vicinity of Tc is 
in agreement with previous observations for other HTSC 
in the low-fluence regim o^°i^^'^'^i^^ and it has been in- 
terpreted as the manifestation of a ex 1/A divergence 
predicted by several theoretical calculations for BCS 
superconductorsiSi^i^^. Experimentally, this divergence 
is smeared out because of the finite superconducting tran- 
sition amplitude. However, from the analysis of the SCS 
decay time increase just below Tc, we can extract quan- 
titative information about the superconducting gap. 

A very useful model to interpret the non-equilibrium 
dynamics of superconductors in the low-intensity regime 
is the Rothwarf-Taylor (RT) modeli^^^s^ (See Ap- 
pendix). In this model, two QPs recombine to form a 
Cooper pair emitting a boson with energy larger then 2A. 
As the reverse process is also allowed, the QP and the bo- 
son populations are in quasi-equilibrium and the actual 
relaxation is determined by inelastic processes. Within 
this model, one can write down a set of coupled rate 
equations, which have analytic solutions in two very im- 
portant limiting cases, the weak and the strong bottle- 
neck regimeai^i^. In the first case, the boson inelastic 
decay rate is fast and the relaxation dynamics is equiva- 
lent to simple bi-molecular dynamics^i^^. In the second 
regime, namely the strong bottleneck one, the inelastic 
decay of the boson population strongly slows down the 
relaxation process. For a given superconductor, the dy- 



namical regime is determined by the particular type of 
bosons considered in this dynamics. However, in both 
regimes and far enough from the critical temperature^, 
a very simple formula for the QP decay rate 7 is valid 



7(T, A{nph),nph) = {uph + nT)T{T, A{nph)) (!) 

where ut are the thermal QPs, Uph are the photo- 
injected ones and r(T, A) is a function of the microscopic 
probabilities for the scattering events involving QPs and 
boson&i^. Given the QP population densities, Uph and 
ut, one can extract from the experimental QP decay 
rate 7 the T(T, A) function in the low-excitation limit. 
We stress that the use of this formula does not imply 
any assumption on the particular boson involved in the 
non-equilibrium dynamics. 

We use the well-known result obtained by Kabanov et 
al.^- that, for a d-wave superconductor with a 2D Fermi 
surface with nodes and for temperatures fc^T << A ~ 
bkBTc (Ref. [33), the QP population at thermal equilib- 
rium, riT, has the form 



riT = lMN{Q){kBTf/A 



(2) 



where N(0) is the density of states at the Fermi level 
and A is the gap value at equilibrium. The temperature 
dependence of ut obtained from Eq. ([2]) is plotted in 
Fig. [5^. To validate Eq. ^, we estimate ut from the 
temperature dependence of the B amplitude (inset of Fig. 
[1]) through the formulaii^ 



B(T) 



2Nph + Uph 



^ 



\'6n\ 



SriT 



(3) 



where Nph is the photoexcited boson population den- 
sity and we assume, similarly to Refs. [13,26, the total 
population density {uph + 2Nph) to be constant in tem- 
perature since the pump fluence is constant. The B{T) 
amplitude in the low-temperature limit is measured at T 
= 10 K. Good agreement is found between ht predicted 
by Eq. ^ and the value obtained from the experimental 
data through Eq. (O (Fig. [2}i). Thus in the following 
we will use the value of nr calculated through Eq. ([2]). 

In Fig. [2^ the experimental SCS decay rate 7 and the 
total QP density, Ntot—nph+nT , are compared. The con- 
tribution due to Uph is estimated in the low-temperature 
limit ■^'^. 

In the high-tcmpcraturc limit we observe that: i) the 
thermal population ht is dominating on Uph', ii) "fexp is 
not following the Ntot temperature dependence. This 
finding suggests that the decay time increase observed 
close to Tc is related to a decrease of r(T, A). We set 
the separation between the high- and low-temperature 
regimes at 60 K, i.e, when ht > riph- Our conclusions are 
independent of the particular choice of this temperature. 

We now verify that the increase in decay time when 
approaching Tc is related to a real divergence arising 
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In this region, we notice a power-law dependence, 
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with an exponent /? ^ 1/2 (solid line in Fig. [5]d), 
which is the same mean-field critical exponent expected 
for the order parameter A in a BCS superconductor. 
In a d-wave superconductor with Tc of 78 K, the su- 
perconducting gap dependence is well approximated by 
A oc (1 — T jTcY^"^ 1 in the temperature range from 60 
to 78 K (see the inset in Fig. [2b). This assumption 
is still a good approximation in the case of overdoped 
HTSC— £"— . A deviation from this exponent is found a 
few degrees below Tc, since Eq. ([T|) is not applicable in 
the close vicinity of Tc^. 

Within the approximate analytic solution (Eq. ([1]))^, 
we can easily derive the power-law dependence of the 
r(r. A) function in the RT approach : 



r(r. A) (X (1 - TjTcT oc A^" 



(6) 



where a '^ 1 with an estimated uncertainty of about 
30%. This power- law dependence can hardly be derived 
by first-principles, particularly if the nature of boson in- 
volved in the pairing mechanism is unknown. 



FIG. 2; Panel a) shows the excitation density for both thermal 
and photoinduced QP, respectively m and n^h (solid lines), 
and their sum Ntot (dashed line) as a function of temperature. 
For the thermal QP we assumed Eq. [2] and a d-wave like 
A(r) dependence and we compared the predicted value with 
the one obtained through Eq. [3] from the experimental B{T) 
amplitude (inset of Fig[l|. On the right axis we report the 
experimental decay rate obtained from the fit in Fig. [l] The 
b) panel shows 7/T''^ as a function of the reduced temperature 
in a double- logarithmic plot. The long-dashed line represents 



the (1 - TjTc 



0/2 



power-law dependence. In the inset we 



show the result of the numerical integration of the d-wave 
BCS gap equation as a function of temperature compared to 
the (X — TjTcY dependence^. In all the panels the vertical 
short-dashed divide the low from the high temperature regime 
(see text). In graph b) the error bars are within the black 
circles size. 



from the fact that r(r. A) — > when A ^^ and we find 
the power-law that controls this divergence. In Fig. [^b, 
we plot the quantity ^exp/T"^ as a function of the reduced 
temperature (distance from the critical temperature) on 
a double- logarithmic scale. Using Eqs. ([1]) and ([5]), we 
find that this quantity is proportional to 



1_ 

rp2 



F(r,A) 



(4) 



in the temperature region where nx is the dominant 
term in the QP density, i.e., above 60 K. 



B. Discontinuity in the fluence dependence 

In Fig. 13^ andlSb: we report the AR/R traces at 10 K 
obtained at increasing pump intensity. Similarly to Sec. 
nil A[ all the curves were fitted using two exponential 
functions convoluted with the time shape of the laser 
pulse. The results of the fit are superimposed to the 
experimental curves. Fig. [3}; shows the SCS amplitude 
and decay time for each fit as a function of pump fluence. 

The low-excitation regime, i.e., the regime where the 
SCS is proportional to the pump fluence is reported in 
Fig. H^ and in Fig. ^ below 26 /xJ/cm^ = $c- In 
agreement with previous works-^— 12^, we assume that 
the SCS is proportional to the photoinduced QP density, 
Tiph- Thus we conclude that Uph increases linearly with 
the pump fluence in this regime. In the zero-fluence limit, 
where Uph — >■ 0, we observe a divergence of the decay 
time (Fig. |3t). This is related to Eq. ([1]), as the total 
QP density, (nph+nr), becomes extremely small at low 
fluence and at low temperature (10 K). 

Above ^c the SCS has a sub- linear dependence and we 
identify this regime as the high-excitation regime^i^ (Fig. 
Wp and Fig. [3t for $ > $c)- In this regime, the SCS 
exhibits a saturation at a critical fluence, ^cr, of ~55 
/iJ/cm^. This saturation means that no more Cooper 
pairs can be destroyed above '^cr and it is considered as 
the evidence of the superconducting condensate vapor- 
ization during the laser pump pulse^i^. The fact that the 
crossover between the linear and the saturated regime 



does not show an abrupt discontinuity here can be justi- 
fied by the spatial profiles of the pump and probe light 
pulses^. The occurence of a real PIPT in this regime has 
been proved by measuring the emergence of a fast com- 
ponent above ^cr in underdoped Bi2212 single crystals^. 
Similarly we observe in Fig. [3b an enhancement of the 
positive signal A associated with the fast e-ph free carri- 
ers relaxation. 

We can compare the SCS in the case of the photoin- 
duced (Fig. [3fc) and the thermally induced phase transi- 
tion (inset of Fig. [1}. In the former case, the parameter 
setting the level of perturbation of the system is the pump 
fluence, while in the latter case this role is played by the 
sample temperature. Within this analogy, the critical 
fluence ^cr, at which the SCS exhibits the saturation, is 
the counterpart of critical temperature, Tc. 

According to non-equilibrium superconductivity 
modelsii^i^ a superconducting-to-normal-state PIPT 
can be either a first-order non-thermal phase transi- 
tion (/i* model)^'^ or a quasi-thermal second order 
phase transition (T* model)^i^. In the T* model the 
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FIG. 4: We report the |A_R/_R| traces in logarithmic scale at 
several pump fluences at T = 10 K, where the exponential A 
component has been subtracted. The dashed lines represent 
the exponential fit of the decay. The solid line refers to the 
curve at <& = 77/iJ/cm^. 
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FIG. 3: AR/R signal, solid lines in panel a) and b), at 10 K 
at different pump fluences. The fit to the experimental curves 
are shown with dashed lines. The fit parameters (decay time 
tb and amplitude B) are reported in panel c) as a function of 
pump fluence. In panel a) the low-intensity regime is shown, 
corresponding to "I> < ^c in panel c). The short-dashed line 
in panel c) is the linear flt in the low-fluence regime. The 
panel b) and the panel c) for "1> > ^c correspond to the high- 
excitation regime. The long-dashed line in panel c) represent 
the saturation value. The error bars In panel c) are within 
the circles size. 



quasi-thermal condition is applied to electrons and high- 
frequency bosons, both being distributed with statistics 
at an effective temperature T*. In both cases the satu- 
ration of the SCS reflects the impossibility of exceeding 
a critical QP density in the stable superconducting 
phasei"— i^. 

In analogy to the thermal case (studied in Sec. 
nil Ap . the quasi-thermal photoinduced vaporization 
would cause the superconducting order parameter to van- 
ish at ^cr- Thus a diverging decay time should be mea- 
sured, according to Eq. (jH]) and Eq. ((H), and in analogy 
to experimental results about the PIPT in CDW— i^^ and 
SDW— compounds. On the contrary, the experimental 
results here reported show a decay time that remains fi- 
nite and below 3.2 ps, thus excluding the quasi-thermal 
origin of the PIPT and suggesting a finite gap at the 
threshold fluence. 

However, the connection between the decay time of the 
SCS and the superconducting gap is firmly established 
only in the low-fiuence regime within the analytical re- 
suhs of the the RT model^^ (gee. ^UTM Eq. ©). To 
extend this concept to the high-fluence regime, when the 
non-equilibrium superconducting gap A{t) could strongly 
vary in time, a validation of the RT model is necessary. 

Before we discuss this point in Sec. IIII C| we report 
two important experimental facts about the dynamics of 
the SCS (See Figg]), 

a) the relaxation dynamics is well reproduced by expo- 
nential decays (dashed lines) in the high-fluence regime, 
as in the low-fluence one. This decay is compatible with 
the RT model, where the gap is assumed to be constant 
in time, thus suggesting that the variations of A{t) are 
small 



b) for $ > $cr, the decay dynamics collapse into a sin- 
gle curve, indicating that the non-equilibrium gap A(t) 
reaches its minimum value at ^cr-, remaining the same at 
higher fiuences. 



C. Rothwarf- Taylor model in the high fluence 
regime: a time-dependent approach 

When a superconducting system is strongly perturbed 
through an ultra-short laser pulse, we expect the super- 
conducting order parameter A to have strong variations 
in time. The parameters /3, 77 and jesc of the RT model 
(see Appendix for definitions) are affected by these vari- 
ations (Refs. Il0ll43( ) and can vary in time. Therefore the 
high-perturbation limit requires a new time-dependent 
Rothwarf- Taylor model. Herewith below we use the fol- 
lowing assumptions, 

i) the time-dependent non-equilibrium superconduct- 
ing gap A{n{t)) can be expressed as a function of n{t) 
considering the T* and /i* model3ii>^>^. In both cases the 
normalized A(n(i)) depends on (l-an(t)) (being n{t) the 
QP density and a a conversion factor) for an s-wave gap 
symmetry and (l-an(t)'^") for a d-wave gap symmetry. 

ii) /?(A(n(i))) is constant, being its temperature de- 
pendence very weak, as reported on YBCO^^ 

iii) r]{A(n(t))) is determined with a fit of the rise time 
at low fluence (See Appendix) and is set to a constant. 
This is a justified by the weak temperature dependence 
of the pair-breaking timei^'l- observed in pump-probe ex- 
periments at low fluencei^ii^ (see Sec. IIII A|) 

iv) 7esc(A(n(t)) is the only time-dependent parameter 
and is responsible for the gap-dependence of T(T, A), ev- 
idenced in Sec. IIII Al fEg. [5]). The value corresponding 
to the unpertubed gap, 7esc(0), can be determined by 
a fit of the experimental decay at low fluence using the 
time-indepedent RT model. (See Appendix) 

v) following a well-established trend in the 
literature^i^ii^ii^ we assume that the cuprates are 
in the strong-bottleneck regime. We can thus conclude 
that r ex 7esc (Ref. SI). 

Under these approximations, the time-dependent RT 
equations are: 



with 



Iphit) 



n = lQp{t) + 2'qp - Pt? 
rip + l3n^/2 - jesc{t) ■ (p - Pt) 



lescit) = JesciO) ■ (A(n(i))/ A(0))2" 

where for an s-wave gap symmetry we have: 

lesc{t) = 7e.c(0) • (1 - an(i))2" 

and for a d-wave gap symmetry: 



(7) 
(8) 



(9) 



(10) 
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FIG. 5: \AR/R\ of the SCS at $ = 26 ^ J/cm^ and the corre- 
sponding fit (solid line) obtained with the numerical solutions 
of the time-dependent Rothwarf- Taylor equations (Eqs. (O 
and ((8])) considering a d-wave gap symmetry. The a opti- 
mal value is (1.0±0.2)xl0~^° cm^. For comparison we also 
show the fits when the parameter a is fixed to the values a 
= 2.45x10"^" cm^ (short-dashed) and a = (long-dashed), 
which is equivalent to the time-independent RT model. In the 

inset we show the corresponding function (7esc(i))^° , which 
is proportional to the non-equilibrium superconducting gap 
A(t) (Eq. |9]). We normalized the gap to the low-fluence value 
(See appendix). 



lescit) = 7e.c(0) • (1 - an{tf/^f" 



(11) 



with a being a conversion factor which set the per- 
turbation of the non-equilibrium gap by QPs and a the 
exponent appearing in Eq. [5] determined in Sec. IIII Al fa 
^ 1). The time-dependent RT equations (Eqs. ([7]) and 
([5])) are then integrated numerically and used to fit the 
experimental curves of the SCS. 

In Fig. [S] we report the results relevant to a sin- 
gle experimental curve at $ = $c and T = 10 K, 
where we assume a d-wave gap symmetry (Eq. (fTTj) ). 
The determined a parameter is (1.0±0.2)xl0^^° cm'^. 
The agreement with the experimental data is very good. 
The predicted normalized non-equilibrium gap, A(t) ex 

(7esc(i))^"j shows a minimum of 60% at the maximum 
SCS and an almost complete gap relaxation (85%) after 
a time delay of 5 ps. Similar results are obtained consid- 
ering an s-wave gap symmetry (Eq. (ITU1) V 

For comparison, we report also two fits where a is 
held constant with values a = 0, equivalent to the time- 
independent RT model, and a = 2.45x10"^" cm^, where 
the gap is suppressed by w 95%. In the latter case the 
predicted relaxation dynamics is extremely slow (of the 
order of several tens of picoseconds) and strongly non- 
exponential. Both these characteristics are not compat- 
ible with the experimental curve reported in Fig. ??, 
thus indicating that, below $(7, the non-equilibrium gap 
is never close to zero. 



Above ^c, the universal and exponential decay dy- 
namics mentioned in Sec. IIIIBI suggests a finite non-zero 
gap after 1-2 ps. It is possibile to estimate a lower bound 
for this finite non-equilibrium gap, considering the op- 
timal a parameter and extracting A{n{t))min from the 
maximum SCS for each curve in Fig. S) The result is 
that, at any time instant t where the RT model is appli- 
cable, the condition A(t) > ^A(O) is realized. 



D. Scenarios for a superconducting-to- normal 
dynamical phase transition 

Although the time-dependent RT model discussed in 
the previous section provides a remarkably good descrip- 
tion of the experimental data, we observe that the time- 
independent RT model is also a valid approximation, pro- 
vided an average gap value (corresponding to the value 
predicted by the time-dependent approach at ~ 2 ps 
delay) is assumed. In particular we demonstrated that 
in our experimental conditions we have A(nph, lOK) > 
^A(0, lOK), thus we can estimate the maximum density 

I 



of photoexcited QP, Uph < nT=7qji « v/P (according 
to the inset of Fig. [2b and Ref. [21). This condition in- 
sures the validity of the analytical formulas obtained in 

Sec. HIS 

We recall Eq. ([T]) and, plotting the decay rate as a 
function of pump fluence (Fig. |B^), we observe the ex- 
pected linear dependence at low fluence. The intercept 
of 0.092 ps~^ is related to the finite value of n^- at T = 
10 K. This linear dependence on the pump fluence fol- 
lows the trend observed on underdoped HTSG^— , but 
was not observed clearly on overdoped samples^i^. We 
define an effective photoinduced decay rate, jph, 



jphiT, A{nph),nph) = "/{T, A{nph),nph) 
-7T(T,A(0),np,=0) 



(12) 



that represents the decay rate exclusively due to the 
photoinjected QPs, where the thermal contribution to 
the decay rate, 7t, (intercept in Fig. [B^) has been sub- 
stracted. We neglect the sample temperature dependence 
since in our pump-probe experiment this parameter re- 
mains constant. We thus obtain the formula 



-fph{A{nph) , riph) = nphT{A{nph)) + nT[r{A{nph)) ~ r(A(0))] « nphT{A{nph)) 

I 



(13) 



In the right-hand side of Eq. [T3] we neglected the sec- 
ond term of the sum. This approximation is valid in the 
high-excitation regime when Uph is much larger than ht, 
as long as F is not going to zero, which is suggested by 
the absence of a decay time divergence. 

Using Eq. ^ and the proportionality between Uph and 
AR/R we can rewrite Eq. (TTS)) as 



Iph 
from which we obtain: 



AR 



A 



2a 



(14) 



Iph 
AR/R 



ex A 



(15) 



1 
We thus define the quantity 6 = ( -^^fj^ ) ° , propor- 
tional to the non-equilibrium superconducting gap, A, 
and we assume a ~ 1 (Sec. IIII Ap . 

In Fig. [61d we report the value of this quantity, nor- 
malized to the zero fluence limit value, obtained from the 
experimental curves of Fig. [3] as a function of the SCS 
amplitude. The gap value decreases with AR/R, hence it 
decreases with increasing Uph, as expected^'^'^. On the 
other side, above the experimental threshold, where Uph 
is constant, the non-equilibrium gap remains constant. 



reflecting the correctness of neglecting the T dependence 
in A{nph,T) and jph{T,A,nph). 

The experimental data are here compared to the pre- 
dictions of the /i* model, where a flrst-order phase 
transition is expectedi"— . Within this model the non- 
equilibrium QP population is described by a Fermi- 
Dirac distribution with an effective chemical potential 
/!*, while its temperature remains at the equilibrium 
value T. In Fig. [B}) we plot the analytical results for 
A(n, T = 0)/A(n = 0, T = 0) in the low-excitation hmit 
as a function of the QP density n for both s-wave and 
d-wave gap symmetrjsS,. The dependence in the d-wave 
case is ex n^''^, while it is linear in the s-wave case, as 
reported in the previous section. 

The fi* model predicts a superconducting-to-normal- 
state phase transition at a flnite gap value, when the 
superconducting-state free energy increases above the 
normal-state value. At 10 K this phase transition should 
take place at a relative gap value of «60% and «65% in 
the s-wave and d-wave case, respectively^, (FiglSlD) 

The main features evidenced by the experimental data 
reported in Fig. [5] are: i) a flnite value of Alnmax,T = 
10 K) of about 1/2 of the equilibrium value at the fluence 
threshold i^cr, h) the linear dependence of A{nph,T = 
10 K) with AR/R and thus with Uph- 

The result i) self-consistently conflrms our initial 
assumption of a non-vanishing F function and of 
A{nph, lOK) > 5A(0, lOK). The experimental value ob- 



8 



tained in this work is remarkably close to the fj,* model 
predictions, even if a quantitative comparison would re- 
quire a detailed analysis of the photoinduced QP distri- 
bution and accurate numerical calculations in the high- 
fluence regime^'"^. 

The result ii) suggests an s-wave like dependence of 
the non-equilibrium gap. Even if this finding seems 
in contradiction with the equilibrium gap d-wave sym- 
metry reported for HTSG^, one should bear in mind 
that the non-equilibrium population photoinduced by a 
100 fs laser pulse in a d-wave superconductor is mostly 
peaked in the anti-nodal region^ because of energy- 
momentum conservation constraint o^^i^^ . This leads to 
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FIG. 6: a) Experimental decay rate extracted by the fit of 
the AR/R traces at different pump fiuences (see Fig. [3]). The 
dashed line is the linear fit obtained in the low-intensity limit. 
The error bars are within the circles size, b) We report 5 as a 
function of AR/R. We normafized the value to the zero flu- 
ence limit, which correspond to the equilibrium gap A(T, 0). 
The error bars account both for the 30% uncertainty in the 
determination of a and for the difference between the results 
obtained within the analytic and the numerical solution of 
the RT equations (see text). On the same graph (right and 
upper axis) we show the analytical results of the fi* model 
for the superconductinggap, A, as a function of QP density 
n as reported in Ref. '2 for both s-wave and d-wave cases. 
The threshold values are calculated numerically within the 
fj,* model^ at a finite temperature corresponding to 10 K in 
our experiment. 



a non-thermal effective distribution that can be repro- 
duced by a Fermi-Dirac statistic with anisotropic effec- 
tive chemical potential. 

The results of such anistropic model are equivalent 
to those obtained for an s-wave superconductor. An s- 
wave like gap symmetry was previously used to explain 
the temperature dependence of Uph in the low-excitation 
regime of HTSC^° and led to a strong controversy in the 
pump-probe experiments interpretation^iiS. We remark 
that this result does not imply an s-wave gap symmetry 
in HTSC at equilibrium and it is rather related to the 
excitation process. 

Both these results suggest a strongly non-thermal 
QP population distribution, thus opening two possible 
physical scenarios behind the superconducting-to-normal 
PIPT: i) the first-order dynamical phase transition pre- 
dicted by the fi* model^"— and supported by the time- 
dependent RT model (Sec. IIIIC[) . ii) a second-order tran- 
sition with a quick (< 1-2 ps) gap recovery, where the 
fundamental assumptions of the RT model break down 
on the short timescales in the high-excitation regime. 

The second scenario goes beyond the treatment pro- 
posed in Sec. IIII CI as the phonon bottleneck would 
be avoided on the short timescales due to the highly 
non-thermal QP distribution. The break down of the 
RT equations implies the decoupling of the relaxation 
dynamics of QP and high-frequency bosons, assessing 
the non-thermal character of the PIPT, in contrast with 
the quasi-thermal PIPT predicted by the T* model^i^. 
Within the first picosecond the superconducting state is 
recovered with a finite non-equilibrium gap, as demon- 
strated in this section. 

However the time-dependent RT model reproduces 
very precisely the experimental dynamics at $£" (Fig. 
[SI Sec. IIIICp . thus suggesting the vahdity of this model 
at excitation levels quite close to the PIPT. Thus, while 
the non-thermal second-order phase transition implies an 
abrupt break down of this model between $c and ^cr, 
the non-thermal first-order transition is compatible with 
the time-dependent RT model also above ^cr- 

The results reported in this manuscript are obtained 
assuming the validity of Eq. ([1]), which is an approx- 
imated analytic solution of the time-independent RT 
equations^. We repeated the same procedure starting 
from the numerical solution of the time-indepedent RT 
equations, which is valid in the whole temperature and 
fluence range. However the flnal results of our work are 
unaffected. The error bars in Fig. [5Jd account for the 
deviation of the analytic results from the numerical one 
and for a 30% uncertainty in the determination of a. 

Both the hypothesis proposed in this section have 
a non-thermal character, at variance with the quasi- 
thermal PIPT reported on other systems, such as 
CDWi^ii^ and SDWii. Since the photoexcitation pro- 
cess is the same, this remarkable difference calls for a 
deeper understanding of the non-equilibrium QP and 
bosons thermalization processes in HTSC in respect to 
CDW and SDW compounds. 



IV. CONCLUSION 

We reported pump-probe experiments on an overdoped 
Y-Bi2212 sample at 800 nm. We explored both the low- 
and the high-excitation regimes to study the origin of 
the recently discovered photoinduced vaporization of the 
superconducting condensate2.i§.. We first verified the ex- 
ponent of the power-law divergence at low-excitation in 
the vicinity of Tc- 

In the high-excitation regime, we showed the lack of the 
quasiparticles decay time increase at the fiuence required 
to photoinduce the transition, in contrast to previously 
reported PIPTs on CDWi^^i^ and SDWl"^ systems. At 
that fiuence ($cr ~ 55 /zJ/cm^) we estimated a finite non- 
equilibrium gap value (2 ps after the interaction with the 
pump pulse) of about 1/2 of the equilibrium gap. 

This finding proves the non-thermal character of the 
superconducting-to-normal-state PIPT, leading to two 
possible phase transition mechanisms: i) a first-order 
dynamical phase transition predicted by the /i* non- 
equilibrium superconductivity modeli"— , and supported 
by numerical solutions of an extended time-dependent 
RT model, ii) a second-order transition with a quick gap 
recovery, where the fundamental assumptions of the RT 
model abruptly break down in the high-excitation regime 
on the short timescales (t < 1-2 ps). 

These findings tackle the fundamental question on the 
interaction of an infrared coherent pulse with the super- 
conducting condensate in HTSC at high excitation pump 
intensity. This is a landmark for the growing field of 
pump-probe techniques on HTSG^'^'^'^i"— i^Sii^. Our ex- 
periment defines the maximum pump fiuence {^cr ~ 55 
/iJ/cm^) which still allows to probe the superconduct- 
ing state of Y-Bi2212; it also demostrates that recent 
pump-probe experiments perfomed on Bi2212 at higher 
pump fiuences^"— are dealing with a dynamical compet- 
ing admixture of superconducting, normal and possibly 
pseudogap phases. 

The considerable difference for the PIPT in HTSC as 
compared to CDW^^'^® and SDW^^ opens an interesting 
perspective for a deeper understanding of the fundamen- 
tal distinctions in the non-equilibrium properties of these 
compounds. 
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Appendix: Rothwarf- Taylor equations 

The non-equilibrium dynamics in superconductors 
is usually successfully interpreted^^' ^'^ within the phe- 
nomenological frame of the Rothwarf- Taylor equationa^^. 



P = Iph{t) -riP + Pn^f^ - lesc{p 



'Pt) 



(A.l) 
(A.2) 



describing the density of excitations n coupled to 
phonons, p being the gap-energy phonon density. The 
non-equilibrium QP and phonons are photo-injected in 
the system through the Iqpit) and Iph(t) terms. A 
gaussian temporal profile of Igp (t) and Iph (t) , with the 
same time-width as the laser pulse, is assumed. The 
coupling of the electronic and phonon population is ob- 
tained through a) the annihilation of a Cooper pair via 
gap phonon absorption (prj term) and b) the emission 
of gap phonons during the two-body direct recombina- 
tion of excitations to form a Cooper pair (/3n^ term). 
In the phonon bottleneck regime {ri>"fesc) the excitation 
relaxation is ultimately regulated by the escape rate of 
the non-equilibrium gap-energy phonons {jescip — Pt) 
term, px being the thermal phonon density). The jesc 
value is determined both by the escape rate of the non- 
equilibrium phonons from the probed region and by the 
energy relaxation through inelastic scattering with the 
thermal phonons. The AR/R superconducting signal 
(SCS) is assumed to be proportional to the solution n(t) 
of Eq. (IA.1[) , in agreement with previous works^"— i2£. At 
low fiuence ($ < ^c), the SCS reported in Fig. |4]is satis- 
factorily reproduced by considering only the Iph (t) term, 
i.e. assuming that the pump energy is mainly absorbed 
through excitation of the phonon population. This re- 
sult is in agreement with both theoretical predictions 
within the T* model^ and experimental observations on 
YBCO^i^ and MgBa^. In the fitting procedure we 
assume (3 = O.lcm^/s, as reported in the literatures^. 
The determined free parameters at low fiuence ($ = 
1 AiJ/cm2) are 7 = (4.5±0.5) ps^^ and 7esc=(3.3±G.l) 
ps~^. These value are compatible with both the results 
obtained on LSCO^ and the theoretical estimations of 
anharmonic processes in YBCOi^.. 
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